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The excellent sensitivity and high resolution capability of wide FoV ground-based 
imaging atmospheric Cerenkov telescopes allow us for the first time to resolve the 
morphological structures of pulsar wind nebulae (PWN) which are older and more 
extended than the Crab Nebula. VHE 7-ray observations of such extended nebulae 
(with field strengths below ~ 20 ^G) probe the electron component corresponding to 
the unseen extreme ultraviolet (EUV) synchrotron component, which measures elec- 
tron injection from earlier evolutionary epochs. VHE observations of PWN therefore 
introduce a new window on PWN research. This review paper" also identifies con- 
ditions for maximal VHE visbility of PWN. Regarding pulsar pulsed emission, it is 
becoming clear that the threshold energies of current telescopes are not sufficient to 
probe the pulsed 7-ray component from canonical pulsars. Theoretical estimates of 
pulsed 7-ray emission from millisecond pulsars (MSPs) seem to converge and it be- 
comes clear that such detections with current generation telescopes will not be 
possible, unless the geometry is favourable. 

"To appear in the proceedings of Cherenkov 2005, Palaiseau, France 



1 The definition of pulsar wind nebulae (PWN) 

The term "pulsar wind nebulae" (PWN, or "plerions" - a term coined by Weiler & 
Panagia Q) signifies a 

(1) filled centre or blob-like form; 

(2) a flat radio spectrum of the form a '--^ to -0.3; 

(3) a well-organised internal magnetic field; and 

(4) a high integrated linear polarisation at high radio frequencies. 

Since this original classification, significant progress with PWN research has 
been made, especially in X-rays, so that we can add the following list of typical 
characteristics: 

(5) A torus and jet near the pulsar, with the direction of the jet reflecting the 
direction of the pulsar spin axis and the torus showing an underluminous region inside 
a characteristic scale radius ~ lO^^cm to ^ lO^^cm, believed to be the pulsar wind 
shock radius ||2J; 
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(6) Evidence for reacceleration of particles somewhere between the pulsar light 
cylinder and r^, leading to a hard X-ray spectrum with a photon index ~ 1.5 to 2.0 
nearr^ (for example: - 1.8 for Crab - |3 1, 1.50 for Vela - |4| and - 2.0for3C58- 

0); 

(7) Evidence for synchrotron cooling (spectral steepening) at r > r^, with the 
size of the PWN decreasing towards increasing energies, as seen from the Crab and 
several other PWN. The photon indices of the cooled spectra range between 2.0 and 
2.5; 

(8) Evidence for a correlation between the spindown power and the average spec- 
tral indices of the pulsar and its PWN JS). 

It is now important to understand how the general PWN properties relate to their 
VHE 7-ray properties. We will highlight a few aspects, following a brief discussion 
of the Crab Nebula: 

The Crab Nebula is considered to be the most important prototype PWN, serving 
as a calibration source for hard X-rays to TeV 7-rays and was also the first source 
to be seen at TeV energies Q. The emission of this source is understood to be due 
to synchrotron-self Compton (SSC), with multiwavelength modelling by |81 and {9\. 
The problem with Crab-like PWN from a ground-based perspective is that the strong 
magnetic field strength in such young PWN results in rapid cooling, so that almost all 
energy is radiated as synchrotron emission, leaving a relatively weak SSC component 
[9 1. The only reason why Crab is still bright at TeV energies, is because the spindown 
flux E/cP ^ 10'^* ergs/s/kpc^ is unusually large compared to all other known pulsars. 
The X-ray photons serve as markers for VHE 7-ray emitting electrons. For the lower 
E PWN, the field is much lower so that the electrons radiating extreme ultraviolet 
(EUV) synchrotron photons also inverse Compton (IC) scatter CMBR into the VHE 
7-ray range. 



2 Optimal young VHE 7-ray PWN 

Rapid SNR/PWN expansion within the first 1,000 years would be a key condition 
to reduce the magnetic field strength fast enough so that synchrotron losses are min- 
imised during the early phase of injection while the spindown power is still quite high. 
If To — Po/2Po is the characteristic age at birth (with Pq the birth period and Pq the 
corresponding period derivative), the spindown power at present (for a pulsar braking 
index of n = 3) would reduce with time as 

The age of the pulsar is then evaluated in terms of the ratio t/ro: If the pulsar is still 
relatively young (i.e. t ~ tq), then E is not much less than Eq, whereas the spindown 
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power for old pulsars (t » tq) decreases with time as (t/ro)^. The evolution of 
the 7-ray flux with time is then evaluated by considering a combination of the time 
evolution of E, and the synchroton lifetime of VHE 7-ray emitting electrons: If E'toV 
is the energy of a 7-ray (in units of TeV) as a result of IC scattering on the CMBR, its 
synchrotron Ufetime in a field of strength B^^ (in units of lO^'"^ G) would be 

T, = i.8{kyy)BzlE:,^f. (2) 

However, the synchrotron Ufetime of an electron emitting keV synchrotron photons 
would be correspondingly shorter: 

Tx - 1.2 {kyr)BZl/^E^l{;. (3) 

In regions of our galaxy where the energy density of dust with temperature 25 K 
exceed the CMBR density significantly, lower energy electrons would produce the 
same TeV 7-ray energies, so that the synchrotron lifetime would increase relative to 
the number given in Eqn. (2). Eqn. (2) however breaks down in the extreme Klein- 
Nishina limit. 

Clearly, if the PWN expanded fast enough up to the present time t, with the 
nebular field strength small enough so that VHE emitting electrons can survive the 
early phases of expansion (i.e. > t), we should then have an optimal VHE emitter. 
In this case we can tap most of the energy (in electrons) ejected by the spinning dipole. 
The SSC process would then be inactive (because of the relatively low synchrotron 
brightness), resulting in IC scattering on the CMBR and possibly far-infrared photons 
from galactic dust grains as the dominant VHE 7-ray production mechanism. 

Gaensler et al. [lOJ found evidence for rapid expansion in G0.9H-0.1 and added 
this source to the PWN associated with PSR B 1509-58, forming a class of rapidly 
expanding PWN. These were also the next two PWN which were detected at VHE 
energies LI liil2J . and the former detection was made despite its large distance to the 
galactic centre. Du Plessis et al. [131 was however the first to predict that rapidly 
expanding young PWN (using PSR B 1509-58 as the prototype) should be relatively 
bright at very high 7-ray energies: The latter authors calculated the field strength 
(~ 8 fiG) associated with such rapid expansion in the PWN of PSR B 1509-58 and 
predicted a flux at 1 TeV, which is relatively close to the recently detected 1 TeV flux, 
although the spectral shape is wrong as a result of the poorly constrained X-ray data 
available at that time. 

A recent review on young PWN with outstanding problems is given by Chevalier 
1141 . Of interest in his discussion is the addition of 3C58: For both the PWN of 
PSR B 1509-58 and 3C 58, the total energy in particles is expected to exceed that of 
the magnetic field. Thus, the magnetisation parameter dtot (defined as the ratio of 
magnetic energy density to particle energy density) is much less than unity, in which 
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case IC radiation competes more favourably with synchrotron emission. Note that the 
Kennel & Coroniti 1 15| model does predict a particle dominated wind at the pulsar 
wind shock radius for Crab (cts ^ 0.003 at rg ^0.1 pc), but equiparition is reached 
further downstream near the bright part of the torus, which is not conducive towards 
the survival of VHE emitting electrons. 

The ideal VHE PWN emitters therefore share at least one or more of the follow- 
ing conditions: 

(a) The overall (total) wind magnetization parameter of the PWN ctot should be 
much less than unity (i.e. a particle dominated wind); 

(b) Rapid expansion of the PWN during its early phases of high power input 
from the pulsar results in the survival of the majority of VHE emitting electrons since 
early epochs. This condition, and the former, are formally consistent with synchrotron 
losses being much less important relative to IC; 

(c) The ideal condition (which includes the first two conditions) is to have the 
lifetime of VHE radiating particles comparable to, or longer than, the age of the sys- 
tem (Tj ^ t), surviving even the earlier epochs when the field was stronger, so that 
the total amount of energy in electrons in the PWN is a significant fraction of the max- 
imal rotational kinetic energy of the neutron star 1(1^/2 at birth. Only adiabatic losses 
are then the main source of losses. In this case we do not expect to see an energy de- 
pendence of the PWN size with changing 7-ray energy - a well known phenomenon 
for PWN where the lifetime of particles exceed the age of the system. 

(d) The spindown power at birth should be much larger than the present spindown 
power, so that, assuming the abovementioned conditions hold, relic VHE emitting 
electrons stored in the reservoir since birth may still contribute to the present VHE 7- 
ray flux. In this case we see the integrated spindown power over a time scale into 
the past, which collects electrons from the epoch of much higher spindown power. 

3 PWN visibility: the 'N-zone' empirical approach 

A simplified approach to estimate the 7-ray visibility of PWN is to take the total X- 
ray flux of the PWN, invert the X-ray spectrum to get the electron spectrum, while 
assuming a single constant field strength. The IC flux (as a result of scattering on the 
CMBR and diffuse galactic photon fields) is then calculated. This phenomological 
approach is probably too simplistic for PWN showing significant gradients in the 
magnetic energy density across the emission volume. 

De Jager 1.161 introduced the following quasi-phenomological model-based esti- 
mate of the expected 7-ray flux from G 0.9H-0. 1 , using what we define here as the 'N- 
zone' model approach: Assume the X-ray spectral index Ti and flux normalisation Fi 
are known for N radial shells (determined from X-ray observations), corresponding 
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to radii Ri, i ~ I, . . . , N. Since the confining pressure eventually tends to decelerate 
the flow, the flow speed Vi should reduce with increasing radius. The total number 
of electrons in the last ('N-th') shell should then also be large compared to the to- 
tal, since the electron residence time, cx I/Vat, is maximal for the last shell, leading 
to this outer volume (between radii Rn-i and R^) contributing mostly to the VHE 
emission via IC scattering. Furthermore, even though we are dealing with an un- 
known magnetic field strength profile with Ri, we may assume that the field strength 
does not change inside the last shell, so that we may apply the constant field strength 
('one-zone model') approach to the bulk of PWN electrons, resident in the last shell. 
This should lead to a conservative estimate of the VHE luminosity. One can then 
either use the full numerical method (using the full cross section for IC scattering) 
or the 'pocket calculator' approximations of |16| to 'invert' the X-ray spectrum of 
the 'N-th' zone to give the approximate IC spectrum corresponding to this important 
VHE emitting volume. 

However, care should be taken with this procedure: do the electrons contribut- 
ing to the X-ray domain also contribute to the VHE domain? To identify the same 
population of electrons which produces synchrotron and IC 7-rays via scattering on 
the CMBR, we can use the expression rederived by 1 17J, giving the synchrotron pho- 
ton energy in terms of the IC scattered 7-ray energy i^xeV in units of TeV for a field 
strength of i3_5 in units of 10 /iG: 

E,yn = 70B_5StoV eV. (4) 

It is thus clear that VHE emission from PWN with these typical field strengths corre- 
spond to synchrotron emission in the Extreme Ultraviolet to soft X-ray domain. For 
the IC scattering of 25 K dust photons in regions of the galaxy where this component 
is important |'18' "13" '16], the constant will be nearly ten times smaller, but extreme 
Klein-Nishina corrections will have to be added. 

Wide FoV optical observations of PWN (if possible) are also essential to com- 
plement VHE observations, since the electrons associated with the optical have even 
longer lifetimes than the VHE 7-ray emitting population and probe the time-integrated 
history even closer to the birth of the PWN. GLAST observations of PWN in the 10 
to 100 GeV domain should also be helpful to provide complementary information. 

Finally, one should be careful when extrapolating X-ray spectra into the optical 
domain (even for the 'N-th' shell) since evolutionary effects may create spectral states, 
which are quite different from our extrapolations. 

4 G 0.9+0.1 and 3C 58: A comparative study. 

HESS recently identified the galactic centre region composite SNR G0.9H-0.1 as one 
of the faintest VHE 7-ray sources LI IJ . For this unresolved VHE source (diameter: 
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2 arcmin in radio and X-rays), 1161 defined the outer 'region 3' of the XMM-Newton 
source II19I as the 'N-th' zone and found that a field strength between 10 and 14 fiG 
can explain the VHE spectrum if the X-ray spectrum (Fx ^ 2.3) extrapolates unbro- 
ken into the EUV range, so that we observe a similar photon index in 7-rays. 

Whereas 'region 1 ' contains the much harder uncooled spectral component (close 
to the injection region at rg), 'region 3' in X-rays and the HESS spectrum are consis- 
tent with the synchrotron-cooled equivalent. The fact that we do not see the cooling 
break in the HESS range should place important constraints on the evolutionary his- 
tory of G 0.9H-0. 1 . The reader is also referred to 1 16 1 for a discussion on the lack of a 
spectral break in the HESS spectrum 1 12 1 of the PWN in PSR B 1509-58. 

The northern hemisphere 3C58 is an important member of this class of PWN: 
whereas it shares 'Crab-like' features, its extended nature should place it in the cate- 
gory of G 0.9H-0. 1 and the PWN of PSR B 1509-58, making it an important object for 
VHE 7-ray studies. It is probably related to the Medieval SNR of 1 181 AD, resulting 
in 1201 deriving an initial spin period of ^ 60 ms, given its present period of 65 ms. 

Comparing 3C58 with G0.9H-0.1 yields the following: 

(a) The inferred spindown power of G0.9H-0.1 1.5 x 10'^^ ergs/s - |21 1) is 
comparable to that of 3C58 (2.7 x lO^^ ergs/s - |20|); (b) The 2-10 keV X-ray lu- 
minosities are 5.4 x lO^^dlg ergs/s for G0.9H-0.1 L19J and ~ 1.3 x lO^^'rfla ergs/s 
(derived from 1221 1 for 3C58; (c) The dimensions of 3C58 are 10 x 6 pc-^, which is 
larger compared to G0.9H-0.1, which has dimensions of ^ 5 x 5 pc^; (d) both show 
evidence of synchrotron cooling, with the X-ray photon index increasing towards the 
edge of the PWN; (e) both show an approximate constant X-ray flux in equal radial 
intervals dR, which would hint at similar radial magnetic field strength profiles; (f) 
for both we find that the size of the PWN at the softest X-ray energies is comparable 
to the size of the radio PWN, which would indicate a cooling timescale comparable 
to, or, longer than the age of the PWN/SNR. In the case of the Crab Nebula, this con- 
dition only holds for the radio to far infrared emitting electrons, which do not radiate 
VHE 7-rays. 

Inspecting Figure 2 of f22l for 3C 58 shows that the X -ray and radio sizes are 
equal for photon energies below 1 keV, whereas the X-ray size shrinks with increas- 
ing energy above 1 keV. This means that electrons radiating 1 keV synchrotron pho- 
tons at the radio boundary should have electron lifetimes equal to the 820 year age 
of the SNR, and from Eqn. (3) we see that this corresponds to a field strength of 
Bjq ~ 13 /iG, which is comparable to the field strengths derived for G0.9H-0.1 and 
PSR B 1509-58 from HESS observations rmT TBl. The X-ray flux of 3C58 is 
twice that of G0.9H-0.1, but the major uncertainty in making VHE 7-ray flux pre- 
dictions is the spectral state of the VHE emitting electrons, which would correspond 
to the unseen EUV synchrotron component. One point of concern is the fact that 
the spindown power did not change significantly since birth (assuming the historical 



6 



association with the Medieval SNR is correct), in which case it was not possible to 
collect significant amounts of relic electrons from an earlier epoch of much larger 
spindown power, despite the favourable condition of rapid expansion and negligible 
synchrotron losses. Thus, whereas conditions (a), (b) and (c) appears to be ideal for 
3C 58, condition (d) may prove to be problematic from a detection perspective. 



5 Middle-aged PWN: new progress 

Vela X, the PWN of the Vela SNR, was the catalyst for the concept of PWN evolu- 
tion 1 23 24 1, as well as the interpretation of middle-aged PWN that are offset from 
their parent pulsars |25|: The basic idea is that a SNR shell usually expands into an 
inhomogeneous ISM, resulting in an asymmetric reverse shock crushing back into the 
PWN. The reverse shock originating from the section where the ISM pressure was 
largest, will crush first into the PWN, creating an offset nebula, as seen from Vela X 
1251 . If the electron lifetime is longer than the crushing timescale, we can still see 
relic electrons radiating for relatively long times in a lower B environment. The ra- 
dio and VHE 7-ray emitting electrons will create the appearance of extended shapes 
that are larger than the X-ray shapes, unless the X-ray emitting electrons also have 
relatively long lifetimes. 

ROSAT and XMM observations of PSR B 1823-13 resulted in the detection of 
the PWN G 18.0-0.7, which was found to be offset by a few arcmin south of its par- 
ent pulsar. This prompted Gaensler et al. f2b] to invoke the offset explanation for 
Vela X (as proposed by |25|) as applicable to G 18.0-0.7. This association was fur- 
ther boosted by HESS observations of this region, showing the detection of highly 
significant resolved VHE 7-ray emission, which is also displaced south of the pulsar 
|27|. Whereas background noise in the XMM image makes it difficult to determine 
the exact X-ray size (> 5 arcmin), the VHE 7-ray source (HESS J 1825- 137) size is 
clearly larger (about 0.5 degrees.) The symmetric compact nebular spectrum with a 
photon index of '--^ 1.6 eventually cools towards the asymmetric offset nebula, with 
photon index ~ 2.3. The latter is also consistent with the VHE 7-ray photon index of 
^ 2.4 (interpreted as IC scattering on the CMBR), which hints at the importance of 
evolutionary effects. 

A systematic investigation of middle-aged pulsars at the edges of resolved centre- 
filled VHE 7-ray sources, combined with follow-up radio and X-ray imaging infor- 
mation, should result in the identification of a new class of PWN, with Ground-Based 
Gamma-Ray Astronomy taking the lead in this new direction. GLAST operations at 
its highest energies (where the angular resolution is best) is also expected to make a 
contribution to this field. 
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6 Pulsars: pulsed VHE 7-ray emission 



The possibility of detecting canonical 7-ray pulsars with present 3"^ generation ground- 
based 7-ray telescopes was discussed some time ago (see e.g. 1281 '). but the fact that 
the spectral cutoffs are below the thresholds of the current generation of telescopes 
forces us to wait for GLAST and 4*'' generation ground-based telescopes (hopefully 
operating at thresholds below 30 GeV) to study the physics associated with spectral 
cutoffs in high-B (canonical) pulsars. 

Studies regarding visibility of VHE 7-rays from millisecond pulsars (MSPs) were 
boosted by the pioneering work of Usov 1 29 1 who showed that 7-rays with energies of 
100 GeV may be expected to escape from the magnetosphere of an MSP. A number 
of studies furthermore centre on 7-ray MPS visibility, e.g. |30 31 1, and specifically 
within a GR framework, e.g. I32ll33ll34ll35l (See Table^for a comparison of different 
authors' converging results for curvature radiation (CR) cutoff energies). 

Two factors determine the visibility of pulsars: spectral cutoff energy and ob- 
servable flux. If the cutoff energy is > the threshold energy of the telescope, and the 
flux sensitivity is < Q.QSE/d^, with d the distance, VHE 7-rays from a pulsar should 
be visible. One of the main uncertainties when modelling pulsars, however, is the 
geometry (i.e. magnetic inclination angle x between the magnetic and spin axes, and 
the angle Q of the observer's line of sight with respect to the spin axis) of the pulsar, 
which significantly influences predictions of cutoff energy. 

Historically it therefore appeared as if MSPs would provide suitable laboratories 
for current (then future) 3"^^^ generation telescopes. However, new calculations and 
final calibration of 3'^'^ generation telescopes indicate that 4**^ generation Cerenkov 
telescopes (see e.g. |36|) and GLAST |37| are generally needed in order to observe 
7-ray MSPs, due to the relatively low cutoff energies I32ll33ll38l . 



Reference 


PSR J 0437-4715 


PSR B 1821-24 




xn 


egg'' (GeV) 


Screened 


xn 


e^l,^ (GeV) 


Screened 


m 


all 


< 17 


no 








m 


20; 35 


- 1 - 10 


no 








f39| 


all 


- 1 - 20 


no 








1391 


10 


< 10 


no 








1331 


10 


- 5 


no 


50 


- 43 


yes 


L35I 


20; 35 


- 10 


no 








1401 


35 


- 10 


no 


50 


> 100 


no 


this paper 


10; 20 


< 10 


no 


50 


< 70 


yes 


this paper 


35 


< 1 


no 


all 


- 0.1 - 150 


yes 



Table 1 : Comparing different authors' converging results for peak CR energy, for similar pulsar parameters. 
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6.1 The unscreened case 



The results in this paragraph build on those described in 1341 1391 . (See e.g. 1331 
^5 40 1 for similar studies). We use the general relativistic B- and E-fields in the 
frame corotating with the pulsar 1 4T\ 1421 1441 and only consider the dominant CR 
component of 7-radiation. 




20 40 60 80 20 40 60 80 

Magnetic Inclination Angle ;x (degrees) Magnetic Inclination Angle x (degrees) 



Figure 1: Contour plot of logj^g of the CR cutoff energy in units of eV vs. observer angle and mag- 
netic inclination angle x- The left hand side is for PSR J()437-4715, representing the "unscreened case", 
while the right hand side is for PSR B 1821-24, representing the "screened case". We used pulsar radius 
R = 10^ cm, moment of inertia I = 3 X 10*^ g.cm-^, and mass M = 1.58Af0. 

In ^441 a critical pulsar spindown power inbreak is given. For pulsars with E = 
ind < £^brcak (with fl = 2tt/P, P the period, and / the moment of inertia), no 
screening due to magnetic pair production will take place above the polar cap (PC) 
(within the framework of the Harding-Muslimov PC-type models). We modelled the 
MSP PSR J 0437-47 15 and used the "unscreened" version of the E-field applicable to 
this case. 

Because this model takes the observer and magnetic inclination angles into ac- 
count, it is possible to find the CR cutoff energy as a function of x and The result 
of this calculation for PSR J 0437-4715 is shown in the left panel of Figure^ The 
CR cutoff energies were estimated from i?^((iA^/(i£')-spectra which we calculated 
for different combinations of x and ( in the range [10,85]. 

Maximum CR cutoff energies < 10 GeV (for the parameters listed in the caption 
of figure are obtained when x ^ C- This should be expected, since the observer 
sweeping close to the magnetic axis will sample radiation due to particle acceleration 
by the highest values of the electric potential. Equality does not hold exactly, as these 
maxima are below the y = x fine. This is due to the asymmetric character of the 
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E-field with respect to the magnetic azimuthal angle cj). The energy maxima occur- 
ing at small values of x also agrees with the fact that the E-field has a cos^ x-term 
dominating at low altitudes, with a sin^ x-term coming into play at larger altitudes. 

6.2 The screened case 

We chose PSR B 1821-24 as an example of an MSP with a screened electric potential. 
We used the approximation of |45| and chose the height of the pair formation front 
(PFF) h = Rpc - (17i?Vc)i/2. 

The CR cutoff energy < 150 GeV as a function of x and C for PSR B 1821-24 
is also shown in Figure [2 (right panel). This case seems to be more symmetric with 
respect to the y — x line and reaches the highest energy values at small x due to the 
fact that the E-field oc cos x- 

What is very important from Figure 1 (for both the "screened" and "unscreened" 
cases), is that a MSP is not expected to be visible for current 3"^ generation telescopes 
(within the GR field definition of |21EIEiy-4)), unless the observer sweeps through 
the magnetic axis, while having the angle between this magnetic axis and spin axis 
substantially less than ^ 45 degrees. 

Finally, the detection of "pair-starved" (low spindown) pulsars is very important 
to probe the naked acceleration potential without the complicating screening effects 
due to pair production. It is however clear that we will have to wait for GLAST and 
4*'^ generation telescopes to answer some of the most fundamental questions related 
to rotating astrophysical dynamos. 
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